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We present the Fermilab E791 measurement of baryon - antibaryon asymmetries in the production of A°, S, 
and A c in 500 GeV/c 7r~-nucleon. Asymmetries have been measured as a function of if and p\ over the range 
-0.12 < x F < 0.12 and p% < 4 (GeV/c) 2 for hyperons and -0.1 < x F < 0.3 and p 2 T < 8 (GeV/c) 2 for the A c 
baryons. We observe clear evidence of leading particle effects and a basic asymmetry even at xf = 0. These are 
the first high statistics measurements of the asymmetry in both the target and beam fragmentation regions in a 
fixed target experiment. 



Particle - antiparticle asymmetry is an excess 
in the production rate of a particle over its an- 
tiparticle (or vice- versa). It can be quantified by 
means of the asymmetry parameter 



where N (N) is the number of produced particles 
(antiparticles) . 

Measurements of this parameter show leading 
particle effects, which are manifest as an enhance- 
ment in the production rate of particles which 
have one or more valence quarks in common with 
the initial (colliding) hadrons, compared to that 
of their antiparticles which have fewer valence 
quarks in common. 

Other effects, such as associated production of 
meson and baryons can also contribute to a non- 
zero value of the asymmetry parameter. 

Leading particle effects in charm hadron pro- 
duction have been extensively studied in recent 
years from both the experimental [ [j] and theo- 
retical point of view [ |) . The same type of leading 
particle effects are expected to appear in strange 
hadron production. Although previous reports of 
global asymmetries in A , H and Q hadroproduc- 
tion already exist [ ||, there is a lack of a sys- 
tematic study of light hadron production asym- 
metries. 

From a theoretical point of view, models which 
can account for the presence of leading particle ef- 
fects in charm hadron production use some kind 
of non-perturbative mechanism for hadroniza- 
tion, in addition to the perturbative production 



of charm quarks [ |^] . 

Given E791's 7r~ beam incident on nucleon tar- 
gets, strong differences are expected in the asym- 
metry in both the xp < and xf > regions. 
In particular, as A (or A c ) baryons are double 
leading in the xf < region while both A (A c ) 
and A (A c ) arc leading in the xf > region, a 
growing asymmetry with \xf\ is expected in the 
negative x f region and no asymmetry is expected 
in the positive xf region. S - baryons are lead- 
ing in both the positive and negative xf regions, 
whereas S + are not. Thus a growing asymmetry 
with \xp\ is expected in this case, are both 
non leading, so no asymmetry is expected at all. 

Experiment E791 recorded data from 500 
GeV/c 7r~ interactions in five thin foils (one plat- 
inum and four diamond) separated by gaps of 1.34 
to 1.39 cm. Each foil was approximately 0.4% of 
a pion interaction length thick (0.6 mm for the 
platinum foil and 1.5 mm for the carbon foils). 
A complete description of the E791 spectrometer 
can be found in Ref. [ || . 

An important element of the experiment was its 
extremely fast data acquisition system [ |J which, 
combined with a very open trigger requiring a 
beam particle and a minimum transverse energy 
deposited in the calorimeters, was used to record 
a data sample of 2 x 10 10 interactions. 

The E791 experiment reconstructed more than 
2 x 10 5 charm events and many millions of strange 
baryons. 

Hyperons produced in the carbon targets and 
with decay point downstream the SMD planes 
were kept for further analysis. 



A were selected in the pir~ and c.c. decay 
mode. All combinations of two tracks with an a 
priori Cherenkov probability of being identified 
as a pir~ combination were selected for further 
analysis if the tracks have a distance of closest 
approach less than 0.7 cm from the decay vertex. 
In addition, the invariant mass was required to 
be between 1.101 and 1.127 GeV/c 2 , the ratio of 
the momentum of the proton to that of the pion 
was required to be larger than 2.5 and the recon- 
structed A° decay vertex must be downstream of 
the last target. The impact parameter must be 
less than 0.3 cm if the particle decays within the 
first 20 cm and 0.4 if decaying more than 20 cm 
downstream of the target region. 

S's were selected in the A°ir~ and c.c. decay 
mode and at the same time f2's were selected in 
the A°K~ and c.c. channel. Starting with a A° 
candidate, a third distinct track was added as a 
possible pion or kaon daugthcr. All three tracks 
were required to be only in the drift chamber re- 
gion. Cuts for the daughter A° were the same as 
above, except for that on the impact parameter. 
The invariant mass for the three track combina- 
tion was required to be between 1.290 and 1.350 
GeV/c 2 and 1.642 and 1.702 GeV/c 2 for S and fl 
candidates respectively. In addition, the 3 and fl 
decay vertices were required to be upstream the 
A° decay vertex and downstream the SMD region. 
For fi's, the third track had to have a clear kaon 
signature in the Cherenkov counter. 

From fits to a gaussian and a linear back- 
ground we obtained 2,571,700 ± 3,100 A° and 
1, 669, 000 ± 2,600 A° from approximately 6.5% 
of the total E791 data sample, 996,200 ± 1,900 
H~ and 706, 600 ± 1, 700 3+ and 8, 750 ± 130 ft" 
and 7, 469 ±120 these last four being from the 
total E791 data sample. The final data samples 
for hyperons are shown in Fig. [j]. 

For charm baryons, the five targets were used. 
In most cases, A c 's decayed in air between the 
target foils, and before entering the silicon vertex 
detectors. All combinations of three tracks con- 
sistent with an a priori Cherenkov probability of 
being identified as a pKn and c.c. combination 
were selected for further analysis if the distance 
between A c decay vertex to the primary vertex 
was at least 5 standard deviations, and the in- 
variant mass was between 2.15 and 2.45 GeV/c 2 . 
To further enrich the sample, we required the A c 
to decay at least five standard deviations down- 
stream of the nearest target foil and between one 
and four lifetimes. The A c momentum vector, 
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Figure 1. A° (upper), 3 (middle) and fl (lower) 
invariant mass plots for the final data samples. 
Left side figures are particles, right side are an- 
tiparticlcs. 



reconstructed from its decay products, was re- 
quired to pass within 3er of the primary vertex. 
It was required that primary and secondary ver- 
tex had acceptable \ 2 P er degree of freedom. We 
also required that at least two of the three A c 
decay tracks be inconsistent with coming from 
the primary vertex. The final data sample, fit- 
ted to a gaussian plus a quadratic background 
has 1,025 ± 45 A+ and 794 ± 42 A~. The invari- 
ant mass plot for the pKix combination is shown 
in Fig. |. 

For each baryon - antibaryon pair, an asymme- 
try both as a function of xf and p\ was calcu- 
lated by means of eq. @. Values for N (N) were 
obtained from fits to the corresponding effective 
mass plots for events selected within specific xf 
and p\ ranges. In all cases, well defined particle 
signals were evident. 

Efficiencies and geometrical acceptances were 
estimated using a sample of Monte Carlo (MC) 
events produced with the PYTHIA and JETSET 
event generators These events were projected 
through a detailed simulation of the E791 spec- 
trometer and then reconstructed with the same 
algorithms used for data. In the simulation of 
the detector, special care was taken to represent 
the behaviour of tracks passing through the dead- 
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Figure 2. pK~ir + and pK + ir~ mass distributions 
showing clear A+ and A~ signals in each case 
in < and xf > regions. From fits with 
a Gaussian and a parabolic background we ob- 
tained 122 ± 17 A+ (a) and 92 ± 15 A" (b) in 
the negative x F and 903 ±42 A+ (c) and 702 ±40 
(d) A~ in the positive xp regions. 



selection criteria, the effect of the kaon con- 
tamination in the beam, the contamination of the 
data sample with D and D s mesons decaying in 
the Kmr and KKix modes and the parametriza- 
tion of the background shape. 

Systematic uncertainties are small and negligi- 
ble in comparison with statistical errors for the 
A c asymmetry. However they are not for the hy- 
perons, and are included in the error bars. 
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ened region of the drift chambers near the beam. 
The behavior of the apparatus and details of the 
reconstruction code changed during the data tak- 
ing and long data processing periods, respectively. 
In order to account for these effects, we generated 
the final MC sample into subsets mirroring these 
behaviors and fractional contributions to the final 
data set. Good agreement between MC and data 
samples in a variety of kinematic variables and 
resolutions was achieved. We generated 5 million 
of A , 16.4 million of H, 4.8 million of f2, and 7 
million of A c MC events. 

Sources of systematic uncertainties were 
checked in each case. For hyperons we looked 
for effects coming from changes in the main se- 
lection criteria, minimun transverse energy in the 
calorimeter required in the event trigger, uncer- 
tainties in the relative efficiencies for particle and 
antiparticle, effects of the 2.5% 

K~ contamination in the beam, effects of K° 
contamination in the A° sample, stability of the 
analysis for different regions of the fidutial volume 
and binning effects. For A c 's we checked the effect 
of varying the main 



Figure 3. A°, E and fl asymmetries as a function 
of xf (upper right) and p\ (lower right). The 
asymmetry for xf (jPr) is integrated over all the 
p\ (xf) range of the data set. The left column fig- 
ures show the predictions of PYTHIA/JETSET. 



Asymmetries in the corresponding xf ranges 
integrated over our p\ and in the corresponding 
Pj, range integrated over our xf range are shown 
in Fig. |^ and ^ for hyperons and A c baryons re- 
spectively, in comparison with predictions from 
the default PYTHIA/JETSET. 

We have presented data on hyperon and A c 
production asymmetries in the central region for 
both xp > and xp < 0. The range of xf cov- 
ered allowed the first simultaneous study of the 
hyperon and A c production asymmetry in both 
the negative and positive xf regions in a fixed 
target experiment. Our results show, in all cases, 
a positive asymmetry after acceptance corrections 
over all the kinematical range studied. 







• ♦ 




1 , T T , I T , , , 1 T , T T I , T T , I , , , , 1 , , , , 1 , , , T 1 , T , , 



-0.1 -0,05 0.05 



I \ \ 



Pt 2 (GeV 2 ) 



0.2 0-25 0.3 





• A c 


- 


o A 








1 














i 


1 


i i I i i t i T t i i i 







-0.1 -0.D75 -0.05 -0.025 0.025 D.05 0.075 0.1 



Figure 4. A+/A c asymmetry as a function of 
xf (upper) and p\ (lower). Full lines are the 
prediction of PYTHIA / JETSET. The asymmetry 
for xf (j>t) is integrated over all the (xf) 
range of the data set. The thin horizontal lines 
are for reference only. Experimental points are in 
the center of the corresponding bin. 



Our results are consistent with results obtained 
by previous experiments [ [| . 

Our data shows that leading particle effects 
play an increasingly important role as \xp\ in- 
creases. The non-zero asymmetries measured 
in regions close to x F — suggest that en- 
ergy thresholds for the associated production of 
baryons and mesons play a role in particle an- 
tiparticle asymmetries. 

On the other hand, the similarity in the A and 
A c asymmetries as a function of xf (see Fig. ||) 
suggest that the ud diquark shared between the 
produced A baryons and Nucleons in the target 
should play an important role in the measured 
asymmetry in the xp < region. However, one 
expects the A c asymmetry to grow more slowly 
than the A° asymmetry due to the mass difference 
between the two particles. 

The PYTHIA/ JETSET model describes only 
qualitatively our results, which in turn can be 
better described in terms of a model including 
recombination of valence and sea quarks already 
present in the initial (colliding) hadrons and ef- 
fects due to the energy thresholds for the associ- 



Figure 5. Comparison between the A° and A c 
asymmetries as a function of xf- The horizontal 
error bars indicate the size of the bin in each case. 



ated production of baryons and mesons [ |^] . 
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